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Abstract 

Effects of gravitational softening on the global structure of self-gravitating 
disks in centrifugal equilibrium are examined in relation to hydrodynami- 
cal/gravitational simulations. The one-parameter spline softening proposed 
by Hernquist & Katz is used. 

It is found that if the characteristic size of a disk, r, is comparable to or 
less than the gravitational softening length, e, then the cross section of the 
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simulated disk is significantly larger than that of a no-softening (Newtonian) 
disk with the same mass and angular momentum. 

We furthermore demonstrate that if r < e/2 then the scaling relation 
r oc e 3 / 4 holds for a given mass and specific angular momentum distribution 
with mass. Finally we compare some of the theoretical results obtained in 
this and a previous paper with the results of numerical Tree-SPH simulations 
and find qualitative agreement. 

Key words: galaxies: kinematics and dynamics - galaxies: structure 
- methods: numerical 
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1 Introduction 



In pure N-body as well as gravitational/hydrodynamical, particle-based sim- 
ulations the gravitational field of the individual particles is commonly soft- 
ened, primarily to suppress effects of two-body gravitational interactions - 
see, e.g., Evrard (1988), Hernquist & Katz (1989), Sommer-Larsen, Vedel & 
Hellsten (1996, paper I in the following). 

In paper I we demonstrated that one has to be cautious when comparing 
results of such simulations with reality when gravity has been softened. This 
was illustrated by considering, as an example, non-rotating, self-gravitating, 
isothermal spheres in hydrostatic equillibrium showing that by introducing 
gravitational softening the structure of such isothermal spheres can be dra- 
matically changed relative to the Newtonian case. This occurs, in qualitative 
terms, when the radial extent of the spheres is less than or of the order the 
gravitational softening length, e (in paper I as well as in this paper the one- 
parameter spline softening proposed by Hernquist & Katz (1989) is used). 

It also follows from paper I that for isothermal spheres of a given mass and 
temperature the size increases with increasing e. In particular, for isothermal 
spheres of characteristic scale r < e/2 (and given mass and temperature) it 
follows from paper I that r oc e 3//2 . 

In this paper we study the global structure of self-gravitating disks in cen- 
trifugal equillibrium for softened gravity. We show that if the characteristic 
size r of the disks is less than or of the order e, then the size of the disks is 
significantly larger than for the Newtonian case for a given mass and angular 
momentum. We also show that if the characteristic size of the disk r < e/2 
then r oc e 3 / 4 (for a given mass and specific angular momentum distribution 
with mass). 

In section 2 we describe the effects of gravitational softening on the global 
structure of self-gravitating disks and discuss the results obtained. Finally 
section 3 constitutes the conclusion. 
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2 Effects of gravitational softening on the global 
structure of self-gravitating disks in cen- 
trifugal equilibrium 



The surface density of a stellar, galactic disk is in general well approximated 
by a truncated exponential (e.g. Freeman 1970, Van der Kruit 1987) 



where S is the central surface density, Rd the exponential scale length, R t 
the truncation radius (i? t ~ 4i? rf ) and R the radial coordinate in the plane 
of the disk. Consider first, for a given mass and angular momentum, the 
one-parameter family of self-gravitating, infinitely flat disks in centrifugal 
equillibrium described by eq. [1] (with R t = 4i?^) for various values of the 
gravitational softening length e. As e increases and the gravitational field 
of the disk becomes increasingly softened one would expect Rd to increase 
(for a given mass and angular momentum) and hence the characteristic cross 
section of the disk 



compared to the Newtonian gravity value <r C) o - in eq. [2] < R 2 > denotes the 
mass weighted average of R 2 . In Figure 1 o" c /<t Ci0 is shown as a function of 
the characteristic size of the disk R rms in units of e. As can be seen from the 
figure, when the characteristic size of the disk is equal to the gravitational 
softening length, a c is about five times larger than <t Cj o due to the effect of 
gravitational softening. For R rms — 0.4e this ratio has increased further to 
about two orders of magnitude. In the limit R rms <C e, <j c /<t c ^ oc (R rms /e)~ 6 , 
as can be easily verified. 

Clearly, gravitational/hydrodynamical simulations where self-gravitating 
disks of size comparable to or smaller than the gravitational softening length 
occur have to be interpreted with care. 

In the simple calculations above the mass and angular momentum was as- 
sumed fixed, but clearly the specific angular momentum distribution, j(M), 
where M is the cumulative mass within radius R, changes with e. Indeed, 
for Newtonian gravity and a given self-gravitating disk in centrifugal equi- 
llibrium with surface density S(_R;e=0) and corresponding specific angular 




(1) 



a c = irR 2 , 



where R rms =< R 2 > 1/2 



(2) 
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momentum distribution j(M), a disk solution E(i2; e) with the same mass 
and j(M) does not necessarily exist for softened gravity with a given non-zero 
softening length e. 

For a disk with surface density given by eq. [1] for Newtonian gravity 
(e = 0, E = Rd — 1, Rt — 4:Rd) we were able to obtain disk solutions 
E(i2; e) with the same mass and j(M) for e > Rd using a numerical, iterative 
algorithm. The solutions are shown in Figure 2 as e 3 / 2 E as a function of 
e -3/4^> ^ gee b e i ow ) f or R d / t = 2~ n ,n = 0,2, ..,8. The corresponding values 
of (x c /cr Cj o are plotted in Figure 1 as filled circles for n — 0, 1, .., 7. 

For < e ^ ^ no solutions could be obtained - we comment further on 
this below. 



2.1 The linear approximation and self-similarity 

By an argument similar to the one given in paper I it can be shown that for 
a system of linear size less than about half the gravitational softening length 
the gravitational acceleration at R < e/2 is approximately linear in R: 

9e (R) = - q R + 0{R 3 ), for R < e/2 , (3) 

where 



AGM, 



g = ^' ^ 

G being the gravitational constant and Md the total mass of the disk. In this 
limit it is possible to obtain a disk solution E(i2; e) given a disk E(_R; 0) with 
a corresponding j(M): 
By definition 

j(M) = Uc (i2)i2 , (5) 
where the circular speed at R is given by (Binney & Tremaine 1987) 

f poo roo 1 1/2 

v c (R) = R 1/2 2nG / [/ J (kR')Z(R ';0)R 'dR'] J l (kR)kdk , (6) 
Jo Jo . 

where Jo and J\ (= — J' Q ) are Bessel functions of the first kind, and 

M(R) = 2ir [ R Z(R'; 0)R'dR' . (7) 
Jo 
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For later use 



dj dj dR 1 u dv c (R) 

[v c {R) + R — ] . (8) 



ciM ciRciM 2tt j RS( j R; 0) L v y d# 
Now, in the linear approximation (assuming e 7^ 0) 

= v c (R) = q l ' 2 R , (9) 



so 



and 



Hence 



j{R) = q^R 2 , (10) 
dj dj dR 2q l ' 2 R q 1 ' 2 



dM dRdM 2vr J RS( J R;e) 7rE(R;e) 



(11) 



s(i?; e ) = (^) 1/2 e - 3/2 (^)" 1 , (12) 

where, by eq. [10], 

RU)=q- 1/A 3 1/2 = (^) 1 'V / 'j 1/2 ■ (13) 

Given a E(i2; 0) eqs. [6]-[8], [12] and [13] can be used to determine the 
corresponding S(_R; e). It is easy to show that the solutions are self-similar: 

Z(R;ae) =a- 3 / 2 Z(a- 3/4 R;e) . (14) 

For E(i2; 0) given by eq. [1], e 3/2 S(i?;e) function of e 3 / 4 i? is shown 
in Figure 2 as a solid line. The solutions obtained previously for R^/e — 
2~ n ,n = 0,2, ..,8, as displayed in Figure 2, converge towards this limiting 
solution as Rd/t — > 0, as expected. 

Given that in the linear approximation a solution E(i2; e) with the same 
mass and specific angular momentum distribution as the Newtonian disk 
E(i2; 0) can be obtained (eq. [12]) it is not surprising that such solutions 
in general can be obtained for Rd/e less than some critical value. But, as 
described in the previous subsection, one should not expect to be able to 
obtain such solutions for all values of Rd/e- 

In the simplified Tree-SPH galaxy formation simulations described in 
paper I self-gravitating, quasi-isothermal gas systems of characteristic size 
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r < e/2 occur at t = 3.2 Gyr in three, otherwise identical, simulations with 
e = 1.5, 3.0 and 6.0 kpc respectively. As discussed in paper I the systems 
become increasingly pressure supported relative to rotational support as e 
increases. 

For systems of equal mass and temperature it follows from eqs. [29] and 
[31] of paper I that one would expect r oc e 3//2 for pressure supported systems. 
Likewise it follows from eq. [14] above that for rotationally supported systems 
of equal mass and j(M) one would expect r oc e 3 / 4 . In Figure 3 we have 
plotted the half-mass radii of the most massive, self-gravitating system at t 
= 3.2 Gyr in the three simulations as a function of e. As can be seen from 
this Figure ri/ 2 {t) increases approximately as e 3//4 as e increases from 1.5 to 
3.0 kpc and approximately as e 3 / 2 as e increases further from 3.0 to 6.0 kpc 
in good agreement with the increasing degree of pressure support found in 
paper I. 

3 Conclusion 

In hydrodynamical/gravitational simulations involving gravitational soften- 
ing an adverse effect of the softening is to make simulated, self-gravitating 
disks artificially large by a significant factor if the characteristic size, r, of 
the simulated disks is less than or comparable to the gravitational softening 
length, e. In such cases the outcome of the simulations should be interpreted 
with due care. 

If r < e/2 then the scaling relation r oc e 3 / 4 holds for disks with a given 
mass and specific angular momentum distribution j(M). 

Finally, the theoretical results obtained in this paper and paper I are 
compared with the results of numerical, Tree-SPH simulations and qualitative 
agreement is found. 
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Figure captions 



Figure 1: The ratio of the characteristic cross section of an exponential disk 
for softened gravity to that of a Newtonian exponential disk for a given mass 
and angular momentum (solid line). The dots show the same quantities for 
solutions of given mass and specific angular momentum distribution j(M). 

Figure 2: The surface density distribution in the linear approximation of 
a disk of the same mass and specific angular momentum distribution, j(M), 
as a Newtonian, exponential disk with R t = 4i?^, displayed as e 3 / 2 S vs. 
e -3/4^> (solid ii ne ). Other curves: Solutions of the same mass and j(M) for 
R d /e = 2~ n ,n = 0,2,. .,8 also displayed as e 3/2 S vs. e~ 3 / 4 R. The solutions 
converge monotonically toward the limiting solution shown by the solid line 
as n increases. 

Figure 3: Half-mass radius as a function of e for the most massive, dense 
gas clumps at t — 3.2 Gyr occurring in the simulations described in paper I. 
Also shown are power law relations of logarithmic slope 3/4 and 3/2 respec- 
tively. 
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